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Abstract Bats have increasingly been recognized as the nat¬ 
ural reservoir of severe acute respiratory syndrome (SARS), 
coronavirus, and other coronaviruses found in mammals. 
However, little research has been conducted on bat 
coronaviruses in South Korea. In this study, bat samples 
(332 oral swabs, 245 fecal samples, 38 urine samples, and 
57 bat carcasses) were collected at 33 natural bat habitat sites 
in South Korea. RT-PCR and sequencing were performed for 
specific coronavirus genes to identify the bat coronaviruses in 
different bat samples. Coronavimses were detected in 2.7% 
(18/672) of the samples: 13 oral swabs from one species of the 
family Rhinolophidae, and four fecal samples and one carcass 
(intestine) from three species of the family Vespertiliodae. To 
determine the genetic relationships of the 18 sequences ob¬ 
tained in this study and previously known coronaviruses, the 
nucleotide sequences of a 392-nt region of the RNA- 
dependent RNA polymerase (RdRp) gene were analyzed phy- 
logenetically. Thirteen sequences belonging to SARS-like 
betacoronaviruses showed the highest nucleotide identity 
(97.1-99.7%) with Bat-CoV-JTMC15 reported in China. 
The other five sequences were most similar to MERS-like 
betacoronaviruses. Four nucleotide sequences displayed the 
highest identity (94.1-95.1%) with Bat-CoV-HKU5 from 
Hong Kong. The one sequence from a carcass showed the 
highest nucleotide identity (99%) with Bat-CoV-SC2013 from 
China. These results suggest that careful surveillance of 
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coronaviruses from bats should be continued, because animal 
and human infections may result from the genetic variants 
present in bat coronavirus reservoirs. 
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Introduction 

Coronaviruses are large, pleomorphic, enveloped viruses con¬ 
taining a single linear, positive-sense single-stranded RNA 
molecule. Coronavirus genomes are approximately 27- 
32 kb in length, the largest continuous RNA genomes among 
mammalian viruses [1]. Coronaviruses are classified into four 
genera: alpha, beta, gamma, and delta [2]. Coronavimses are 
the second most prevalent cause of the common cold in 
humans [3]. In livestock and poultry, coronavimses are recog¬ 
nized causes of enteric and respiratory infections that are often 
fatal in young animals [4]. An outbreak of severe acute respi¬ 
ratory syndrome (SARS) in 2002 and 2003 resulted in infec¬ 
tion of 8096 people worldwide, with 774 (9.5%) of them 
dying from this novel human coronavims [5]. SARS corona- 
vims is believed to have been acquired from an animal spe¬ 
cies, and the Chinese horseshoe bat may have been the orig¬ 
inal source of SARS infection [6-9]. The emergent Middle 
East respiratory syndrome coronavims (MERS-CoV) may al¬ 
so have originated in bats, with fragments of MERS-CoV 
genes being identified in bats from both Saudi Arabia [10] 
and Africa [11]. More than 200 vimses have been isolated 
from and detected in bats [12]. The ability of bats to fly and 
migrate, as well as the large sizes of their social groups, pre¬ 
disposes bats to acquire and maintain vimses [13]. 
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Given the import of MERS into South Korea [14] and the 
presence of SARS in the relatively close geographic location 
of China [9] (Fig. 3), together with the fact that bats are a 
reservoir for coronaviruses, the prevalence of coronavirus in¬ 
fection in Korean bat species should provide valuable infor¬ 
mation. Although SARS-like and MERS-like bat 
coronaviruses have already been reported in the feces of 
Korean bats [14], the prevalence and presence of 
coronaviruses in other bat samples have not been described. 
Therefore, we investigated the prevalence of coronaviruses in 
Korean bat species using various kinds of samples, including 
oral swabs, fecal and urine samples, and carcasses. The 
coronaviruses found were characterized for their relationships 
to each other and to coronaviruses isolated in South Korea and 
in other countries. 

Materials and Methods 
Sample Collection 

A total of 672 samples (332 oral swabs, 245 fecal samples, 
38 urine samples, and 57 bat carcasses) were collected at 33 
sites of natural bat habitat from January 2016 to September 
2016 (Table 1). Bats were captured by net for collection of 
oral swabs, fecal samples, and urine samples, and released 
immediately after sampling. Fecal samples were also col¬ 
lected from guano. The bat carcasses collected were natu¬ 
rally dead bats found at the site where the other samples 
were taken. Swab samples were kept in viral transport me¬ 
dium at 4 °C, and organs collected from carcasses and other 
samples were stored at -80 °C before processing. 

Genome Detection and Sequencing 

All fecal and tissue samples were resuspended in 1% 
antibiotic-antimycotic solution (Corning, USA) diluted in 
phosphate-buffered saline (PBS), and clarified by centrifu¬ 
gation at 3500xg for 10 min. RNA from 200 qF sample was 
extracted with the QIAamp® Viral RNA Mini Kit (Qiagen, 
Germany) and eluted in 60 pB RNase-free water. cDNA was 
synthesized by PrimeScript First Strand cDNA Synthesis 
Kit (TaKaRa, Japan) following the manufacturer’s instruc¬ 
tions. Bat-CoV screening was performed by a pan- 
coronavirus PCR method based on primers used by Poon 
et al. (Corona 1 forward, 5'-GGTTGGGACTATCC 
TAAGTGTGA-3', and Corona 2 reverse, 5'-CCAT 
CATCAGATAGAATCATCATA-3') [15], followed by se¬ 
quencing of the amplified product to confirm Bat-CoV iden¬ 
tification. The pan-coronavirus primers were used to ampli¬ 
fy and sequence a 440-bp segment of the highly conserved 
RNA-dependent RNA polymerase (RdRp) gene. Cycle se¬ 
quencing reactions were performed using the BigDye® 


Terminator Cycle Sequencing kit version 1.1 (Applied 
Biosystems, Foster City, CA, USA). Reactions were puri¬ 
fied using SigmaSpin™ Post-reaction Clean-Up Columns 
(Sigma, St. Fouis, MO, USA) and sequenced on an ABI 
PRISM 3130 automated capillary DNA sequencer 
(Applied Biosystems) according to the manufacturer’s in¬ 
structions. Primer sequences were excluded from the raw 
sequence data, and the sequences we analyzed were 392- 
bp long. All sequences obtained in this study were submit¬ 
ted to GenBank (accession numbers KY432454-432471). 

Phylogenetic Analysis 

The nucleotide sequences were aligned and compared to 
52 selected human and animal CoV sequences available 
from the GenBank database using ClustalW software im¬ 
plemented in BioEdit version 7.0.9.0. The phylogenetic 
trees were drawn using the neighbor joining method using 
the maximum composite likelihood model with MEGA 7 
software. 

Results 

Prevalence of Coronavirus Infection in Bats 

Coronaviruses were detected in 13 oral swabs, one carcass 
(intestine), and four fecal samples from a total of 672 sam¬ 
ples, giving an overall detection rate of 2.7% (18/672). 
Four of the 14 bat species tested were found to harbor 
coronavirus. Thirteen positive oral swabs were detected 
from one species of the family Rhinolophidae , which 
accounted for 24.3% (163/672) of the samples, and the 
other positive samples were from Vespertilionidae bats that 
comprised 55.5% (373/672) of the samples. 

The 13 coronavirus-positive oral swabs were collected 
from bats at an abandoned mine in Jeonbuk province (Fig. 
1). Oral swabs and other samples (n = 60) were obtained from 
three species of bats, Rhinolophus ferrumequinum , 
Miniopterus schreibersii , and Myotis macrodactylus , but 
coronaviruses were only detected in samples from 
R. ferrumequinum (Table 2). One carcass was collected from 
a cave in Chungbuk province (Fig. 1). In total, 94 samples 
were obtained from six bat species (M schreibersii , Myotis 
petax , Vespertilio sinensis , Murina leucogaster , 
R. ferrumequinum , and M. macrodactylus ), and coronavirus 
was detected from a V. sinensis sample. Four coronavirus- 
positive fecal samples were collected in four different regions, 
in Gyeongbuk province (Andong, Yeongdeok, and Gyeongju) 
and in the metropolitan city Gwangju (Fig. 1). In the forest of 
Andong, 43 samples in total were collected from three species 
of bats ( Eptesicus serotinus , M. petax , and Pipistrelus 
abramus ), and coronavirus was detected from one 
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Table 1 Bat sample information 

Region 1 

Region 2 

Site name 

Site type 

Major bat species 

No. of samples 









Oral swab 

Feces 

Carcass 

Urine 

Total 

Gangwon 

Yeongwol 

YW 1 

Forest 

Eptesicus serotinus 

11 

18 



29 

Gyeongnam 

Changwon 

CW 1 

Abandoned mine 

Rhinolophus ferrumequinum 


7 



7 

Gyeongbuk 

Gyeongju 

GJ 1 

Abandoned mine 

Rhinolophus ferrumequinum 

70 

2 

6 

12 

90 



GJ 2 

Forest 

Myotis aurascens 

1 

1 



2 



GJ 3 

In area of private houses 

Pipistrelus abramus 


10 



10 




Forest 

Pipistrelus abramus 


5 



5 


Andong 

AD 1 

Forest 

Eptesicus serotinus 

25 

18 



43 



AD 2 

Forest 

Eptesicus serotinus 

6 

2 


2 

10 


Yeongdeok 

YD 1 

Abandoned mine 

Myotis macrodactylus 

2 


3 


5 



YD 2 

Forest 

Pipistrelus abramus 


5 



5 


Yeongju 

YJ 1 

Abandoned mine 

Myotis nattereri 

23 

6 

1 


30 



YJ 2 

Forest, stream 

Pipistrelus abramus 


5 



5 


Yeongcheon 

YC 1 

Abandoned mine 

Miniopterus schreibersii 



12 


12 


Cheongsong 

CS 1 

Forest, valley 

Myotis petax 

4 

1 



5 

Gwangju 

Shinchon 

GwJ 1 

Forest 

Eptesicus serotinus 

11 

6 


1 

18 


Seokgok 

GwJ 2 

Cave 

Rhinolophus ferrumequinum 

5 




5 

Daegu 

Dalseong 

DS 1 

Forest, stream 

Pipistrelus abramus 


5 



5 

Ulsan 

Ulju 

US 1 

Abandoned mine 

Rh inolophus ferrumequinum 

12 


2 

8 

22 



US 2 

In area of private houses 

Eptesicus serotinus 

1 

1 



2 




Forest 

Eptesicus serotinus 

31 

31 



62 



US 3 

Forest 

Myotis aurascens 

9 




9 

Jeonnam 

Gurye 

GR1 

Forest 

Eptesicus serotinus 

12 




12 


Muan 

MA 1 

Abandoned mine 

Myotis macrodactylus 

7 



1 

8 



MAU 

Cave 

Myotis formosus chofukusei 

2 

5 



7 


Jindo 

JD 1 

Abandoned mine 

Rhinolophus ferrumequinum 

10 

1 

1 


12 


Hampyeong 

HP U 

Cave 

Myotis formosus ch ofukusei 

3 




3 

Jeonbuk 

Sunchang 

SC 1 

Abandoned mine 

Rhinolophus ferrumequinum 

43 

3 

10 

4 

60 

Chungnam 

Seosan 

SSU 

Abandoned mine 

Myotis formosus ch ofukusei 

3 

4 



7 

Chungbuk 

Danyang 

DY 1 

Cave 

Miniopterus schreibersii 

7 

65 

22 


94 



DY 2 

Cave 

Myotis macrodactylus 


2 


2 

4 



DY 3 

Cave 

Miniopterus schreibersii 

18 

1 


6 

25 



DY 4 

Cave 

Miniopterus schreibersii 


25 



25 



DY 5 

Forest, valley 

Myotis ikonnikovi 

1 




1 



DY 6 

Forest, valley 

Myotis ikonnikovi 

8 

6 


2 

16 



DY 7 

Forest 

Murina leucogaster 

7 

10 



17 

Total 





332 

245 

57 

38 

672 


P. abramus bat. In Yeongdeok, a total of five samples was 
collected from P abramus living in the forest, and one sample 
was positive for coronavirus. In the forest of Gyeongju, five 
samples in total were collected only from P abramus bats, and 
one of the samples was positive for coronavirus. In the forest 
of Gwangju, 18 samples in total were collected from two 
species of bats (E. serotinus and Myotis aurascens) and coro¬ 
navirus was detected from one E. serotinus bat (Table 2). 


Phylogenetic Characterization of Bat Coronaviruses 
Isolated in Korea 

To determine the genetic relationships between the 18 bat 
coronaviruses obtained in this study and previously known 
coronaviruses, the 392-nt RdRp sequences were analyzed 
phylogenetically. All the sequences obtained in this study 
belonged to the betacoronavirus genus (Fig. 2). The sequences 
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Gangwon 


Seoul 

Gyeonggi 


Chungbuk 


Chungnam 


Daejeon 


Jeonbuk 


Gyeongbuk 

A 


Daegu 


Gwangju 


Gyeongnam 



Ulsan 


Busan 


Jeju 


Fig. 1 Location of South Korea and geographical location of coronavims-positive bat samples. Black circle Bat-CoV from oral swabsm. Black up- 
pointing triangle Bat-CoV from fecal samples. White up-pointing triangle Bat-CoV from a bat carcass 


from 13 oral swabs (16B0122, 129, 132, 133, 134, 135, 139, 
141, 143, 145, 150, 154, and 155) were grouped with the 
SARS-related bat betacoronaviruses, lineage B. Sequences 
from a carcass (16BT3) and four fecal samples (16BF104, 
109, 211, and 244) were grouped with MERS-like bat 
betacoronaviruses, lineage C (Fig. 2). 

The 13 sequences from oral swabs were 97.1-100.0% 
identical with each other, and the closest known CoV strain, 
with 97.1-99.7% nucleotide identity, was Bat-CoV-JTMC15 
that was recently isolated from R. ferrumequinum in China. 
The nucleotide identity with SARS-CoV was only 86.7- 


89.0%. One sequence from a carcass showed the highest iden¬ 
tity (99%) with Bat-CoV-SC2013 that was isolated from 
Vespertilio superans (synonym of V. sinensis ), also from 
China [16]. The nucleotide identity with MERS-CoV was 
83.9%. Four sequences from fecal samples shared 96.9- 
100.0% identity with each other and showed 94.1-95.1% 
identity with Bat-CoV HKU5 [17]. The nucleotide identity 
with MERS-CoV was 84.1-84.4%. 

Due to shorter reported sequences from Kim et al. [14], 
another phylogenetic tree was constructed from the 355-bp 
RdRp sequences of Korean bat coronaviruses obtained by 
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Fig. 2 Phylogenetic tree 
constructed from 392-nt RdRp 
sequences of Korean bat 
coronaviruses detected in this 
study and from other 
coronaviruses. Phylogenetic trees 
were constructed by using the 
neighbor joining method and 
bootstrap values were determined 
by 3000 replicates. Scale bar the 
estimated genetic distance of 
these viruses. Black circle se¬ 
quences from oral swabs Black 
up-pointing triangle sequences 
from fecal samples. White up- 
pointing triangle sequences from 
carcass 
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a-CoV 


0.05 


Kim et al. [14] and in this study to analyze the genetic relation 
between Bat-CoVs in South Korea (Fig. 3). Thirteen se¬ 
quences from oral swabs were clustered with Bat-CoV B15- 
21, which was detected in fecal bat samples collected from an 
abandoned mine in Gangwon province. The major bat species 
present there was R. ferrumequinum , but the coronavirus¬ 
positive bat species were not specified [14]. Nucleotide iden¬ 
tity between Bat-CoV B15-21 and the 13 clustered sequences 
was 97.7-100.0%. The other three Korean bat coronaviruses 
(Bat-CoV B15-8, B15-40, and B15-41) belonged to the 
alphacoronavims genus, but we did not detect coronaviruses 
belonging to this group in our study. 


Discussion 

There are approximately 24 species of bats in South Korea, 
and most are insectivorous and relatively small [18]. Bats in 
South Korea usually do not migrate over distances greater 
than 1000 km [19], but apart from two species found only in 
South Korea, they are generally distributed in neighboring 
countries such as China, East Siberia, Japan, and Taiwan. 
Although it is not clear if bats go between China and other 
countries, some species that migrate over 500 km can carry 
pathogens to South Korea. The fact that the 18 bat 
coronaviruses (Bat-CoVs) detected in this study were all 
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Fig. 3 Phylogenetic tree of Korean bat coronavimses detected by Kim 
et al. [13] and in this study, based on 355-nt sequences. Phylogenetic trees 
were constmcted by using the neighbor joining method and bootstrap 
values were determined by 3000 replicates. Scale bar the estimated 


genetic distance of these viruses. Black circle sequences from oral swabs. 
Black up-pointing triangle sequences from fecal samples. White up- 
pointing triangle sequences from carcass 


closely related to the Bat-CoVs identified from bats in China 
suggests the possibility that bat coronaviruses have been in¬ 
troduced from China and other countries into South Korea, or 
vice versa. This also can indicate the migration of bats over 
borders between countries. 

In this study, the prevalence of coronavims in Korean bat 
samples was 2.7% overall. The prevalence by sample type was 
5.4% (13/332) for oral swabs, 1.8% (1/57) for carcasses, 1.6% 
(4/245) for fecal samples, and 0.0% (0/38) for urine samples. 
Reported rates from other countries for coronavirus prevalence 
in bats were 6.5% (64/985) and 5.3% (50/951) for China, 8.2% 
(16/195) for Italy, and 5.3% (32/606) for Mexico [20-23], not 
greatly different from our figure of 2.7% for prevalence in 
South Korea. 

Fourteen out of 24 Korean bat species were screened for 
Co Vs in this study. The species positive for Bat-CoVs were 
R. ferrumequinum , V. sinensis , P abramus , and E. serotinus. 
The Bat-CoVs that showed the highest nucleotide identity with 
our isolates were obtained from the same species that we iden¬ 
tified as positive, except for E. serotinus. To our knowledge, 
this is the first time a Bat-CoV has been detected from 
E. serotinus in Asia. In Europe (Italy), a lineage C 
betacoronavims was detected from E. serotinus in 2014 [24]. 


Bat-CoVs were usually detected at higher rates in alimentary 
samples (fecal samples, rectal swabs, and intestine) than in 
respiratory samples (oral and throat swabs) [7, 15, 17, 22, 25, 
26]. We collected alimentary and respiratory samples and urine 
samples. However, Bat-CoVs were not detected in urine sam¬ 
ples. Animal coronavimses usually produce either enteric or 
respiratory infections. Animal models show similar clinical fea¬ 
tures to SARS and MERS patients; the majority of whom pre¬ 
sented with respiratory symptoms but in some cases also suf¬ 
fered from enteric complications [27]. Unlike other human co¬ 
ronavims infections, a number of MERS cases were associated 
with acute renal failure. MERS-CoV replication in the human 
kidney suggests the potential of shedding vimses in urine [28]. 
Bat-CoV isolate Bat SL-CoV-WIVl (KF367457) showed in- 
fectivity to Rhinolophus sinicus kidney cell lines, which means 
it may be possible to detect Bat-CoV in urine [6]. 

The current study revealed that Korean bats have SARS- 
and MERS-related Bat-CoVs, consistent with earlier work 
[14]. The previous study detected alpha- and 
betacoronavimses only from fecal samples in four regions, 
two in Gangwon, one in Chungbuk, and one in Gyeongbuk 
[14]. The only betacoronavims RdRp sequence Bat-CoV 
B15-21 (KU528590.1) was clustered with 13 Bat-CoVs 
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(16B0122, 129, 132, 133, 134, 135, 139, 141, 143, 145, 150, 
154, and 155) detected from Danyang, Chungbuk province. 
The identity between the 14 sequences was very high at 97.7- 
100%, even though the sample types and collection sites were 
different. The species of bat infected with Bat-CoV B15-21 
was not specified in the paper, but the major bat species of the 
sample collection site was R. ferrumequinum ; the same spe¬ 
cies in which we detected 13 Bat-CoVs. The same Bat-CoVs 
were detected from fecal samples in the prior study [14] and 
from oral swabs in our study. Spike proteins, which are trans¬ 
lated from spike genes, define viral tropism by receptor spec¬ 
ificity and by membrane fusion activity during entry into cells 
[29]. Although the spike genes were not identified in this 
study, full spike gene analysis of Bat-CoV B15-21 suggested 
a low potential for human emergence. One MERS-like 
betacoronavirus, Bat-CoV B15-1-6 (KU528586.1, 
KU528585.1), was identified by Kim et al. [14], but the sim¬ 
ilarity of this isolate to our isolates could not be analyzed 
because the RdRp sequence was not obtained. High nucleo¬ 
tide identity between Bat-CoV B15-1 -6 and Bat-CoV-SC2013 
(KJ473821) suggests that Bat-CoV B15-1 -6 and 16BT3 might 
be very similar. 

Bat-CoVs have been identified in many countries including 
the Australasian region [10, 11, 14—16, 20-26, 30, 31]. But 
Bat-CoVs identified in Korea were genetically close to Bat- 
CoVs from China, a close neighbor of Korea. According to 
the results of Bat-CoV research in China, the degree of host 
restriction for coronaviruses in bat populations was high. For 
example, similar viruses were detected in the same species of 
bat in different regions, approximately 1600 km apart, while 
two different bat species in the same habitat had different 
coronaviruses [22]. In our results, three Bat-CoVs 
(16BF109, 211, and 244) from P. abramus were clustered 
with each other, even though samples were collected from 
different habitats in different regions (Andong, Yeongdeok, 
and Gyeongju). The Bat-CoV detected from E. serotinus bat 
in the metropolitan city Gwangju (16BF104) was also clus¬ 
tered with these three Bat-CoVs mentioned above. However, 
there was no matching virus for 16BF104 found in the same 
species of bat in other locations, nor in other species in the 
region where samples were collected. Our data could support 
the possibility of relationships between genetic variation, geo¬ 
graphic locations where Bat-CoVs were isolated, and bat host 
species, but more data on Bat-CoVs in Korea would be need¬ 
ed to address this. 

In summary, the prevalence of coronavirus in Korean bats 
was determined and found to be comparable to or lower than 
the rates in other countries. We detected SARS-related and 
MERS-like Bat-CoVs that were genetically very similar to 
Bat-CoVs identified in China. For the study of Bat-CoV sea¬ 
sonal occurrence and transmission, continuous monitoring 
during all seasons is required, and ideally alimentary, respira¬ 
tory, and other samples need to be collected from each bat. 


Future work to determine the complete genome sequences of 
Bat-CoVs from South Korea would give a more complete 
picture of their pathogenicity and the possibility of infectious 
spread to other animals and humans. 
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